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Since the 1950s, cellulose pyrolysis has been the subject of intense study, with kinetic analyses 
forming a major part of these studies. They represent useful tools for a better understanding of 
the physicochemical process and for the proper design of industrial pyrolysis units. Until recently, 
the methods most frequently used in these analyses were based on model-fitting, i.e. the fitting 
of the experimental data to a number of mathematical models. Nowadays, other methods, so- 
called “model-free” methods, are considered to be more suited. These axe based on the principle 
that, at constant conversion, the reaction rate depends only on temperature. In its first part, this 
short review presents the particularities and drawbacks of the traditional model-fitting models. 
Subsequently, several main contributions in this field are listed and discussed. Finally, the more 
suited “model-free” (isoconversional) methods are explained and several main studies presented, as 
well as a comparison of this method with the model-fitting ones. 
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Introduction 

Cellulose is the most abundant biopolymer on 
earth. It is found in the stalks of vegetal species and 
in wood. Cellulose is usually isolated from wood by 
an acidic treatment, affording a wood pulp, of which 
the insoluble fraction is denoted as “a-cellulose”. It is 
the type of cellulose in most frequent use in the phar¬ 
maceutical and food industries. A high content of this 
compound is found in the white cotton fibres (around 
94-96 mass %) which are considered to be pure cel¬ 
lulose. Cotton cellulose is used for manufacturing mi¬ 
crocrystalline cellulose (MCC) and filter paper. 

Together with hemicellulose and lignin, cellulose is 
an important compound of biomass, representing from 
50 mass % up to 75 mass %. This is one of the rea¬ 
sons behind the intensive studies of cellulose as a case 
model for wood and biomass pyrolysis. Nowadays, py¬ 
rolysis is a promising environmentally-friendly tech¬ 
nical process for reduction of greenhouse gases emis¬ 
sions and valorisation of wastes. The amount of CO 2 
released during biomass pyrolysis is equal to that ab¬ 
sorbed by the plant during its growth. Despite these 


advantages, in 2008 there were only 30 to 50 high ca¬ 
pacity commercial pyrolysis plants in existence world¬ 
wide. This stemmed from the drawbacks represented 
in particular by the endothermic character of the pro¬ 
cess and the necessity of treating the solid and liquid 
fractions obtained for a more adequate valorisation. 

An extensive literature exists on cellulose pyrol¬ 
ysis, as pyrolysis represents a promising way of ob¬ 
taining fuels from renewable sources. Thus, in 2010, 
in Germany, energy consumption from biomass repre¬ 
sented 7.9 %, out of a total of 11 % of renewable energy 
(AGEE-Stat, 2011). However, little progress has been 
made up to the present and the real mechanism has 
not yet been found. 

This review summarises the main contributions on 
the subject of cellulose pyrolysis kinetics published in 
the open literature. Two directions are emphasised: 
the model-fitting vs the so-called “model-free” meth¬ 
ods. First, the drawbacks of the classical model-fitting 
techniques are outlined. Next, the main reasons for re¬ 
garding the “model-free” methods as more suitable are 
considered. A comparative discussion of the two meth¬ 
ods is also presented. However, it should be noted that 
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this short review covers partially only the existent lit¬ 
erature on the kinetic analysis of cellulose pyrolysis at 
different heating rates (from slow to flash). 

Problems regarding thermo-analytical 
kinetic modelling of solids decomposition 

At the beginning of the 21st century, Galwey 
(2004) emphasised the major problem with kinetic 
studies based on thermal analysis experiments. 
Among his conclusions, the lack of physical meaning 
for the Arrhenius parameters calculated and the im¬ 
possibility of finding a real kinetic mechanism for the 
process from thermal analysis data alone are worth 
noting. According to the same author, supplementary 
tests (X-ray analysis, microscopy, etc.) are required. In 
the thermo-analytical literature at that time he also 
identifiedthe existence of ambiguities in the definition 
of essential terms (“mechanism”, “rate constant”, “ac¬ 
tivation energy”) and a lack of order in the results pre¬ 
sented. Little progress has yet been made in resolving 
these problems. 

Accordingly, the main drawback of this kind of ap¬ 
proach, dealing with the Arrhenius parameters cal¬ 
culated, remains unsolved. As Vyazovkin and Wight 
(1997) stated, “if we are to interpret them in terms 
of the transition-state theory, they are not applica¬ 
ble to solid state reactions”. In the solid state reac¬ 
tions, due to the stable and tight array-structure of the 
sample, the energy distributions functions of Maxwell- 
Boltzmann are not suited. However, other energy dis¬ 
tribution functions can be used such as the Fermi- 
Dirac function for electrons and the Bose-Einstein 
function for photons. In this case, E a is related to the 
activation enthalpy and A is related to the frequency 
of lattice vibrations, as also stated by Vyazovkin and 
Wight (1997). The physical significance of these pa¬ 
rameters refers to the elemental processes of solids, 
i.e. nucleation, nuclei growth, etc. Meanwhile, from 
thermo-analytical experiments, the general kinetic pa¬ 
rameters are obtained. Another major problem of this 
kind of approach is the empirical nature of the kinetic 
models tested (White et al., 2011). Further investiga¬ 
tions are needed in order to identify the real mech¬ 
anism of solids’ degradation, due to the fact that,in 
this case, the thermal changes (including the chemical 
changes) are often more complex than is recognised. 

All these problems represent a consequence of the 
origins of modern thermal analysis kinetics, as noted 
by Galwey (2004). These origins are located in a spe¬ 
cialised branch of chemistry concerned with the ther¬ 
mal decomposition of solids, known as “crystolysis re¬ 
actions” . The theory applied in these cases was based 
on geometric models applicable to heterogeneous reac¬ 
tions in crystals. For these crystals, the stoichiometries 
were regarded as already well established. 

In conclusion, for solid state thermal degradation 
kinetics, studied using the thermal analysis methods, 


the values of the Arrhenius parameters calculated de¬ 
scribe a given step of the process in a general manner 
and have a different physical meaning than in the case 
of gas or liquid phase transformations. 

Main contributions to cellulose pyrolysis 
model-fitting kinetic models 

Several models describing the cellulose pyrolysis ki¬ 
netics have been published with different degrees of 
accuracy. As Di Blasi (1993) stated, the studies of cel¬ 
lulose pyrolysis largely followed the chemical regime 
conditions (intrinsic kinetics) while several attempted 
to model the process under ablative regime condi¬ 
tions (apparent kinetics). Ablative pyrolysis is con¬ 
tact pyrolysis, allowing large particles to be treated 
as in flash pyrolysis (Lede, 2012). The intrinsic kinet¬ 
ics conditions provide the information for identifying 
the real degradation mechanism of cellulose pyrolysis. 
Despite huge efforts, the real kinetic model remains 
unknown. From the perspective of the chemical trans¬ 
formations, this process is regarded as proceeding in 
three stages, differentiated by the temperature range 
(Garcia-Perez, 2008). Thus, the pyrolysis stage is the 
first thermal degradation observed at temperatures up 
to 700°C. The second stage between 700-850 °C is rep¬ 
resented by the secondary pyrolysis. Finally, tertiary 
reactions generally occur between 850-1000X1. Glob¬ 
ally, the pyrolysis reaction comprises homogeneous 
solid-phase processes, heterogeneous solid-vapour in¬ 
teractions, and also gas-phase transformations. 

Garcia-Perez (2008) proposed a detailed complex 
model for cellulose and biomass pyrolysis, explaining 
the formation of some secondary pyrolysis products 
(PAHs). Furthermore, he identified the main prod¬ 
ucts and their formation mechanisms from each stage. 
Thus, in the primary pyrolysis stage, cellulose-derived 
molecules are formed as fluid species, namely levoglu- 
cosane (LVG) and hydroxylacetaldehyde. From the 
secondary pyrolysis stage, gaseous olefins, aromatics, 
CO, CO2, H2O, H2, and secondary condensed oils 
(phenols, aromatics) are produced. Finally, from the 
tertiary reactions stage, PAHs are formed together 
with non-condensable gases (CO, CO2, H 2 , H 2 0). The 
PAHs formed in the course of cellulose pyrolysis are 
either gases, such as the methyl derivatives of aromat¬ 
ics (toluene, indene, methyl acenaphthylene, methyl 
naphthalene)or condensed gases or tars, like naph¬ 
thalene, acenaphtylene, anthracene/phenantrene, and 
pyrene. Coke and soot are formed in the primary and 
secondary reactions. Evans and Milne (1987) noted 
that coke (“char”) was formed by liquid and organic 
vapours thermolysis, while soot was formed by the 
homogeneous nucleation of high temperature decom¬ 
position products of hydrocarbons in the gas phase 
processes. It is generally accepted that a degradation 
step takes place prior to pyrolysis. Kilzer and Broido 
(1965) were the first to observe this. They explained 
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the evolution of water during this stage as being the 
effect of exothermic cross-linking reactions. The cor¬ 
responding temperature range is approximately 100- 
220 °C, as confirmed in a recent study by Zhu et al. 
(2012) who used an in-situ visualisation technique. 
The same study confirmed that the primary pyrolysis 
step occurred in the temperature range of 230-380 G 
for a final experimental temperature of 500°C. 

On the basis of the above chemical mechanisms 
for cellulose pyrolysis, the existing global kinetic mod¬ 
els for this process were developed. A classification of 
the kinetic models for pyrolysis is given by Di Blasi 
(1993) as follows: (a) one-step global reaction with 
no prediction of product yields dependence on reac¬ 
tion conditions; (b) multi-reaction empirical models 
correlating experimental data with several reactions; 
(c) semi-global models which include primary and sec¬ 
ondary reactions with lumped pyrolysis products (tar, 
gas, and “char”). As noted by Lede (2012), the most 
representative models for cellulose pyrolysis were pub¬ 
lished between 1950 and 1980. The same author states 
that these models are classified into two main cate¬ 
gories: models that take into account the formation of 
an intermediate and models that neglect its formation. 

Fig. 1 presents an historical overview (Capart et 
al., 2004; White et al., 2011; Lede, 2012) of the most 
important model-fitting kinetic models of the cellu¬ 
lose pyrolysis. In this figure, the abbreviations of the 
listed models are given in the grey boxes. The follow¬ 
ing notations were used for the different types of reac¬ 
tions considered in the kinetic models listed, namely 
by their corresponding kinetic constants (see Fig. 1): 


Cellulose —■> Char (1) 

Cellulose —Tar (2) 

Cellulose —> Gases (3) 

Cellulose ———> Intermediate (4) 

Low DP cellulose ———♦ Activecellulose (AC) (5) 
Cellulose —Char + Tar (6) 

Cellulose ——> Char + Gases (7) 

Cellulose —> Char + Tar + Gases (8) 

Intermediate ———> Char (9) 

Intermediate — fcl ° > Tar (10) 

Intermediate ———* Gases (11) 

Intermediate — kl2 > Char + Gases (12) 

Tar ——--- > Gases (13) 

Tar kli > Char + Gases 
Intermediate —- 14 > Char + Gases (BRmodel) 

Tar — —- - > Char + Gases (15) 

Tar fcl6 > Char (16) 

Chari —Char 2 + Gases (17) 


Char 2 —Char 3 + Gases (18) 

A classification of all the models reviewed is given 
below. This considers the linkage or otherwise of the 
final lumped products obtained (names of the mod¬ 
els underlined are those taking into consideration the 
formation of an intermediate, usually AC): 

i) Tars, gases, chars: 

- Model 3R (Agrawal, 1988) with three parallel 
reactions directly from cellulose; 

- Model MKB (Agrawal (1988), “modified Kilzer- 
Broido” model); 

- Model L (Lewellen et al., 1977) is a complex ki¬ 
netic model which explains the decomposition steps 
both in the cellulose matrix and in the ambient atmo¬ 
sphere of the sample; the intermediate is LVG; 

- Model GP (Garcia-Perez, 2008) based on a lit¬ 
erature review regarding the formation of PAHs from 
biomass pyrolysis. 

ii) Char + tar + gases: 

- Model FI (Stamm, 1956) is the first kinetic model 
for cellulose pyrolysis; it is the 1st order kinetic model 
formulated for the pyrolysis data from Douglas fir cel¬ 
lulose; 

- Model VI (Varhegyi et al., 1997) for high- 
pressure conditions and sealed crucibles proving that 
water catalyses the process. 

iii) Char + tar + gases, gases: 

- Model BSZ (Banyasz et al., 2001). 

iv) Char + gases, tar, gases: 

- Models SI (Shafizadeh, 1968) and S2 (Shafizadeh, 
1982) with the difference that SI takes into account 
the transformations of gases into “char + gases” and 
the inverse; 

- Model DKM (Diebold, 1994) considers the sec¬ 
ondary pyrolysis that consists of the degradation of 
tars to gases at high temperatures; 

- Model W (Scott et al., 1989), “Waterloo” pyrol¬ 
ysis, based on flash pyrolysis experiments on Avicel 
cellulose in a fluidised bed reactor operating at 520 G 
and 101 kPa. 

v) Char + gases, tar: 

- Model B (Broido & Nelson, 1975) is the first 
semi-global model; it assumes the formation of an in¬ 
termediate denoted as AC (“active” cellulose or “an- 
hydrocellulose”) and the further degradation of the 
primary “char + gases” into secondary and subse¬ 
quently tertiary “char + gases”; 

- Model KB1 (Kilzer & Broido, 1965); 

- Model BS (Shafizadeh & Bradbury, 1979), “Broi- 
do-Shafizadeh” model, is frequently cited in the lit¬ 
erature for cellulose pyrolysis; it was established from 
pyrolysis experiments under vacuum (2 kPa); this was 
needed to prevent secondary reactions and to improve 
recovery of the primary products; 

- Model V2 (Varhegyi et al., 1997) considers AC 
to contain a mixture of cellobiose, cellotriose, other 
polymers of cellobiose, glucose, and fructose; 




k, (lower DP) ^ volatiles 

*« (tar) k l3 



Cellulose a, solid, ♦(1-a,) gas 

7 k„ a 2 solid 2 + (1-a 2 ) gas 


Cellulose- 1 - Intermediates —♦ Char, H : C, Gas 
, k t k 12 

' -*s-‘Char, Volatiles, H 2 0,Gas 



^ kj -Hydroxyacetaldehyde,Formaldehyde,CO 


- Solid, + vols.—* Solid, + vnls. —- Solid, - vols. 


* Mono- and oligomer type products 


C l I—*HAA, acetol, aldehyde, ketone, acid 

1 J—- Furfural, 5-HMF,etal. 
kit —* Other heterocyclic compound 
—*CO + C0 2 + CH« + C„H m +H 2 
— k n -- Char + H 2 0 + Gas *- 


L \ 


_[__| Jr„ -• Gas, (HOCH 2 CHO,CO,H 2 CO...) 

_____ t L, Intermediate < 

mbs I Celluloses ^ Char, 

[Ti *s H LVG + Cellobiosan I + H 2 0 + C0 2 


0.95 HAA ♦ 0.25 Glyoxal + 0.20 CHjCHO + 
n 0.20 CjH»0 + 0.25 HMFU + 0.20 C0 2 ♦ 

12 10.15 CO + 0.1 CH 4 + 0.9 HjO + 0.65 Char 


H y* eoty^iccM'^et Cellulose < 'Decompositionproducts(glyoxal.HAA...) 

lvw * 7 "char (5 wt.%| ♦ Tar 


- Model LIU (Liu et al., 2008); 

- Model F (Ranzi et al., 2008) is a detailed ki- 
ic model, taking many elementary reactions into 


vi) Char + gases, gases: 

- Model BR (Broido et al., 1973) takes 
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the formation of lower DP cellulose, AC, and LVG as 
intermediates; 

- Model DB (Di Blasi, 1996) introducing a supple¬ 
mentary reaction which takes the secondary pyrolysis 
into account, i.e. the degradation of tars into gases at 
high temperatures. 

vii) Char + gases: 

- Model AF (Alves & Figueiredo, 1989) considers 
further degradation of the primary “char + gases” into 
secondary and subsequently tertiary “char + gases”; 

- Model MBS (Liao et al., 2004) also considers the 
secondary reactions of levoglucosan and gases crack¬ 
ing. 

viii) Char + tar, tar, gases: 

- Model RZ (Ranzi et al., 2008) is a multi-step 
lumped mechanism based on TGA experiments previ¬ 
ously published. 

ix) Char + gases, char, gases: 

- Model H (Hopkins et al., 1984) is based on flash 
pyrolysis experiments. 

x) Char, tars: 

- Model KB (Broido, 1976) is based on experi¬ 
ments realised from TGA under vacuum. 

xi) Char, gases: 

- Model MBY (Mamleev et al., 2007) considers the 
intermediate formation of cellobiosan linked with LVG 
and gases. 

All these models neglect the mass and heat trans¬ 
fer phenomena that become rate-limiting when a high 
mass of the sample and/or high temperatures are used. 
In this case, other models are needed, at intra-particle 
level; this case is not discussed below. Of the models 
listed above, the first order global model (FI) was in¬ 
tensively used for this process in the literature. More¬ 
over, out of the semi-global models, the BS model is 
the most important, and is still considered as a ref¬ 
erence for the kinetic modelling of cellulose pyrolysis, 
albeit contested by Varhegyi et al. (1994). New mod¬ 
els were proposed in the last decade and are consid¬ 
ered better than the BS model or its derivatives, all of 
which have as their main drawback the fact that they 
take into account the formation of an intermediate. 
This was and remains a controversial issue in cellu¬ 
lose pyrolysis, since any such intermediate is difficult 
to identify and quantify. 

The new models consider that a kinetic of nuclei- 
growth best describes the process. The models used for 
describing this phenomenon are either of the Avrami- 
Erofeev or of the Prout-Tompkins type (Capart et 
al., 2004; Mamleev et al., 2009; Barud et al., 2011). It 
is important to note that these nuclei-growth models 
give similar results to the “model-free” methods. In 
this concept, Mamleev et al. (2009) proposed a mech¬ 
anism involving the formation of several melted cav¬ 
ities (nuclei) with an encapsulated liquid inside the 
cellulose’s microfibre. More recently, Sanchez-Jimenez 
et al. (2011) proposed a new model which was better 
than the Avrami-Erofeev, Prout-Tompkins, or other 


nucleation or growth mechanisms. This new model is 
an acceleratory chain-scission model based on break¬ 
age of the polymeric cellulosic chains. In addition, it 
makes no assumptions as to the exact chemical nature 
of the bond breakages. Another type of model applied 
to the kinetic study of this process was a statistical 
model, known as “the distributed activation energy 
model” (DAEM) (Sonobe & Worasuwannarak, 2008). 
This has been in use for biomass kinetic analysis since 
1985 but it has attracted little interest among the sci¬ 
entists from this domain. 

The first order (FI) model is considered by 
Banyasz et al. (2001) as satisfactorily describing the 
process at constant heating rates up to 100 °C min -1 . 
This is due to the negligible effect of the high tem¬ 
peratures pathway under these conditions. Moreover, 
Di Blasi stated that the FI model could be considered 
as rigorously representative of intrinsic kinetics (Di 
Blasi, 1998). However, Volker and Rieckmann (2002) 
observed that the pyrolysis behaviour of small cel¬ 
lulose samples can be described more accurately by 
multi-step rather than by single-step processes. 

Table 1 shows the published values for these pa¬ 
rameters as calculated using the FI model. Those ref¬ 
erences where the values of kinetic parameters var¬ 
ied with the heating rate are omitted from the list, 
since they may be regarded as unreliable on the ba¬ 
sis of the ICTAC 2000 project conclusions that the 
use of isoconversional methods and multiple heating 
rates experiments are valuable tools for describing the 
multi-step thermo-analytical processes (Dickinson & 
Heal, 2009a, 2009b). However, studies of the FI ki¬ 
netic model performed with this type of cellulose are 
few. It can be seen that the use of microcrystalline 
cellulose (MCC) afforded values of these parameters 
which are among the highest (150-240 kJ mol -1 and 
2.5 x 10 17 -5.01 x 10 21 min -1 ). These are superior 
to those obtained using a-cellulose (including filter 
paper) or other types of cellulose. When using a ni¬ 
trogen flow of 60-70 mL min -1 or a helium flow of 
100-150 mL min -1 , it can be observed that the ki¬ 
netic parameters’ values increase with the decrease in 
the final reaction temperature, irrespective of heat¬ 
ing rate, type of cellulose or sample mass. Thus, in 
the case of N2 flow, at temperatures up to 1000 °C, 
the corresponding kinetic parameters’ values were: 
2.82 x 10 14 min -1 and 182.7 kJ mol -1 , while at 
500°C the values were higher: 3.01 x 10 17 min -1 
and 210 kJ mol -1 . For the He flow, at 1000°C these 
values were: 6 x 10 17 min -1 and 210 kJ mol -1 , 
while at 900°C they were higher: 2.88 x 10 18 min -1 
and 226.5 kJ mol -1 . From the studies reviewed us¬ 
ing the FI kinetic model presented in Table 1, it 
can be seen that the calculated energy is situated 
in the range of 82.7-244 kJ mol -1 with a median of 
192.3 kJ mol -1 and the pre-exponential factor range 
is 6.7 x 10 5 -1.28 x 10 20 min -1 with a median of 
5.11 x 10 19 min -1 . 
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Table 2. Published values for kinetic parameters of several main semi-global models of cellulose pyrolysis in inert atmosphere 






/ 


E a 

Reaction 



Reference 

Model 

Year 

— 

— 



constants 

Reaction 

Observations 




n 


- 1 

kj mol - 




Conesa et al. (1995) 

KB1 

1975 

7.9 

X 

10 11 

153.1 

k 2 

Eq. (2) 

35 % char 65 % gas 




1.9 

X 

10 16 

197.9 

fc 7 

Eq. (7) 




1995 

2.3 

X 

10 17 

217.7 







3.2 

X 

10 17 

225.5 




Di Blasi (1998) 

BR 

1998 

4.7 

X 

10 11 

152.8 

fe 4 

Eq. (4) 

Cellulose discs, 0.076 cm thick 




6.6 

X 

10 16 

206.8 

feio 

Eq. (10) 





2.2 

X 

10 17 

210.4 

fcl4 

Eq. (14) 




1979 

1.7 

X 

10 21 

242.4 

&4 

Eq. (4) 

Cellulose powder, 0.5 cm 




1.9 

X 

10 16 

197.9 

kio 

Eq. (10) 

thick; isothermal tube furnace 




7.9 

x 

10 11 

153.1 

ki2 

Eq. (12) 


Conesa et al. (1995) 

BS 

1995 

1.1 

X 

10 24 

274.4 



Filter paper 




6.2 

X 

10 10 

200.9 







4.9 

X 

10 16 

207.5 




Reed and Cowdery (1987) 

DKM 

1987 

1.7 

X 

10 21 

243.0 

fc 4 

Eq. (4) 

_ 




4.0 

X 

10 7 

109.0 

fc 7 

Eq. (7) 





2.1 

X 

10 13 

204.0 

fell 

Eq. (11) 





1.1 

X 

10 5 

61.0 

fcio 

Eq. (10) 





7.8 

x 

10 11 

153.0 

ki2 

Eq. (12) 


Conesa et al. (1995) 

3R 

1988 

3.0 

X 

10 13 

178.9 

ki 

Eq. (1) 

a-Cellulose 




2.0 

X 

10 16 

200.6 

ks 

Eq. (8) 




1995 

1.5 

X 

10 17 

226.9 







3.7 

X 

10 17 

219.2 




Conesa et al. (1995) 

MKB 

1988 

6.4 

X 

10 1B 

195.8 

k 2 

Eq. (2) 

a-Cellulose 




1.8 

X 

10 12 

159.3 

fc 4 

Eq. (4) 





3.4 

X 

10 17 

214.6 

fen 

Eq. (11) 





1.9 

X 

10 16 

203.7 

kg 

Eq. (9) 




1995 

2.5 

X 

10 8 

120.1 







1.6 

x 

10 24 

279.3 







5.4 

X 

10 16 

208.5 







3.8 

X 

10 1B 

206.8 




Conesa et al. (1995) 

AF 

1989 

7.3 

X 

10 13 

174.0 

fc 7 

Eq. (7) 

Filter paper, ai = 0.75 




1.9 

X 

10 14 

178.0 

fel 7 

Eq. (17) 

Filter paper, a 2 = 0.99 




1.2 

X 

10 14 

179.0 

fel8 

Eq. (18) 

Filter paper, <13 = 0.22 



1989 

1.3 

X 

10 11 

139.0 



Cotton wool, ai = 0.87 




4.2 

X 

10 12 

157.0 



Cotton wool, a 2 = 1.00 




1.1 

X 

10 17 

209.0 



Cotton wool, a.3 = 0.16 



1995 

3.3 

X 

10 17 

217.5 



Filter paper, a\ = 0.70 




6.7 

X 

10 11 

128.2 



Filter paper, a 2 = 0.78 




9.7 

X 

10 12 

147.5 



Filter paper, a3 = 0.16 

Miller and Bellan (1997) 

DB 

1994 

1.7 

X 

10 21 

242.4 

fc 4 

Eq. (4) 

_ 




2.0 

X 

10 16 

196.5 

feio 

Eq. (10) 





7.8 

X 

10 11 

150.5 

fel2 

Eq. (12) 





2.6 

x 

10 8 

108.0 

fel3 

Eq. (13) 


Marra (2009) 

KB 

1997 

1.3 

X 

10 16 

167.5 

fc 4 

Eq. (4) 

ni = 0 




1.9 

X 

10 1B 

196.0 

kg 

Eq. (9) 

n 2 = .1.5 




5.6 

x 

10 17 

216.6 

feio 

Eq. (10) 

n 3 = 1.5 

Banyasz et al. (2001) 

BSZ 

2001 

1.0 

X 

10 16 

198.9 

fe 3 

Eq. (3) 

Filterpaper - Whatman 41 




1.4 

X 

10 12 

148.2 

fcs 

Eq. (8) 


Ranzi et al. (2008) 

F 

2008 

4.8 


10 1B 

192.6 

fc 4 

Eq. (4) 

- 
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Table 2 shows the published values for the kinetic 
parameters calculated with several main global multi- 
step models (KB1, BR, BS, DKM, 3R, MKB, AF, DB, 
KB, BSZ, F). A wide scatter of parameter values can 
be observed, confirming the empirical character of the 
approaches used in these kinetic modelling studies. Be¬ 
cause of this variation, no realistic assumptions can 
be made regarding the validity of the kinetic models 
tested. The variation is observed even for the same 
model when tested in different years by the same au¬ 
thor, under the same experimental conditions (KB1, 
BR, 3R, MKB, AF). 

Despite the large number of kinetic studies on cel¬ 
lulose pyrolysis, few insights into understanding of 
the process have been achieved. In general, the pro¬ 
cess chemistry is only globally described by pseudo¬ 
components which consider groups of chemical species. 
In conclusion, these models are important, especially 
for small samples where only intrinsic kinetics occurs. 
In the mid-nineties, Vyazovkin (1996) suggested the 
use of the term “kinetic scheme” instead of “mecha¬ 
nism” in the thermal analysis kinetic studies, as being 
more appropriate. These schemes represent a sequence 
of steps affecting or producing a modification of a mea¬ 
sured physical property. The present author endorses 
this recommendation as suited for further studies in 
this field. 

Main contributions to cellulose pyrolysis 
“model-free” kinetic models 

Since Galwey’s observations (Galwey, 2004) on the 
empirical nature of the kinetic parameters calculated 
by model-fitting methods, more authors have become 
aware of this problem. Thus, the so-called “model- 
free” methods were proposed as a solution. In the last 
decade these new methods have been used more and 
more frequently for the kinetic analysis of solids’ ther¬ 
mal processes. They are based on the principle that, 
at constant conversion, the reaction rate depends only 
on temperature, as can be seen in the mathemati¬ 
cal expression used for their description (Vyazovkin 
& Wight, 1997): 


'din (da/di) 1 _ E a 
d T- 1 \ a _ R 


This is due to the greater sensitivity of the solid phase 
reactions to pressure, temperature, crystal size, etc. 

The isoconversional methods give values with real 
physical meaning for the activation energy; however, 
it is inappropriate to denote them as “model-free”, 
since they only postpone the choice of a certain ki¬ 
netic model until an estimate of the activation en¬ 
ergy has been made, as outlined by Sewry and Brown 
(2002). Moreover, Capart et al. (2004) observed a 
qualitative agreement in the kinetic parameters ob¬ 
tained both from isoconversional and subsequent non¬ 
linear regression methods. Among the authors inter¬ 
ested in this field, Sergey Vyazovkin (Vyazovkin & 
Dollimore, 1996; Vyazovkin & Wight, 1997; Brown 
et al., 2000) and Nicolas Sbirrazzuoli (Sbirrazzuoli et 
al., 2009) made major valuable contributions. They 
outlined in their studies that isoconversional methods 
enable the acquisition of information on the control¬ 
ling mechanisms for a global reaction. Furthermore, 
by applying these new methods, the same mathemat¬ 
ical approach can be used in resolving the system of 
equations. 

The isoconversional methods for the calculation of 
Arrhenius parameters are classified into two major 
types, in respect of the mathematical treatment ap¬ 
plied to the TGA data. Thus, the approaches are ei¬ 
ther differential or integral, based on graphical meth¬ 
ods and they require at least three experiments. 
Among the differential methods, Friedman method is 
the simplest and most frequently used. It is based on 
the rate equation’s differential form (Capart et al., 
2004; Aboulakas & Harfi, 2008): 

ln(,^) = In [/!/(«); (20) 

where fJ/( °C min -1 ) is the heating rate. 

As for the integral methods, the best-known is the 
Flynn-Wall-Ozawa (FWO) method which uses linear 
equations for calculation of E a (Aboulakas & Harfi, 
2008): 


ln/3 
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( 21 ) 


where 


where a is the conversion, t /min is the time, T/K is 
the temperature, E a /(kJ mol -1 ) is the activation en¬ 
ergy, and R /(J mol -1 K -1 ) is the universal gas con¬ 
stant. 

For this reason they are also known as “isocon¬ 
versional” methods. Usually, the E a -a curves enable 
recognition of a multi-step process (Brown et al., 2000; 
Gavillon, 2007; Dickinson & Heal, 2009a, 2009b). This 
can be explained by the fact that, for solid state ki¬ 
netics, unlike liquids and gases, it is possible to have 
a variation of the activation energies with conversion. 




dT = 


/3R P \RT ) 


(22) 


f(a) is the kinetic model of the process and p(E a /RT) 
is Doyle’s approximation for calculating the so-called 
“temperature integral” used in the isothermal model¬ 
fitting methods. 
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However, in this case, the solution of Eq. (21) 
uses several approximations which introduce some er¬ 
rors into the kinetic parameters values. In order to 
overcome this drawback, a non-linear isoconversional 
method was proposed by Vyazovkin and Dollimore 
(1996). This method is based on the minimum con¬ 
dition obtained from the integral form of the general 
Eq. (21): 


v L Ta ’j) fa _ 

/'(^a«;T„ i )A 

* ’ 


(23) 


where I is the identical matrix, n is the total number 
of experimental points and i, j are the numbers of the 
tested heating rates. 

To avoid the use of an approximation for the tem¬ 
perature integral, the method proposed by Budrugeac 
and Segal (2003) can be used. This covers the case 
of solid-gas decomposition using isothermal experi¬ 
ments. 

Other isoconversional methods cited in the litera¬ 
ture are those proposed by Kissinger-Akahira-Sunose 
(Eq. (24)) (Aboulakas & Harfi, 2008) and Li-Tang 
(Eq. (25)) (Li & Tang, 1997, 1999): 
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instruments (Capart et al., 2004). This is the case 
for the FWO (Flynn-Wall-Ozawa) method, the KAS 
(Kissinger-Akahira-Sunose) method and the differ¬ 
ent Vyazovkin methods. Other isoconversional expres¬ 
sions exist in the literature but are rarely used. Ta¬ 
ble 3 gives the literature kinetic parameters obtained 
with these models under different experimental con¬ 
ditions. It may be observed that, irrespective of the 
isoconversional kinetic model used, the heating rates 
(from slow to fast), the type and flow-rate value of 
the carrier gas, and the type of cellulose, the val¬ 
ues are close to each other, especially for the acti¬ 
vation energy. From this type of method, the mean 
value of the activation energy reviewed is very close to 
that obtained from the FI model-fitting kinetic model, 
namely 194.3 kJ mol -1 from the isoconversional mod¬ 
els vs 192.3 kJ mol -1 from the FI model. This demon¬ 
strates that, among the model-fitting kinetic models, 
the FI model is suited to accurately describing the 
kinetics of cellulose pyrolysis. 

The isoconversional methods afford kinetic param¬ 
eter values that can vary with conversion (Vyazovkin, 
1996; Vyazovkin & Dollimore, 1996; Vyazovkin & 
Wight, 1997). In this case, E a represents a temper¬ 
ature coefficient for the global reaction rate, showing 
the sensibility of this reaction rate to temperature. 
The other kinetic parameter, A, represents a scalar 
factor of the global reaction rate and characterises the 
reaction’s amplitude or intensity. 


/ ( ln/ ^) da = / ^ da + G (“> “o) (25) 

where: 

G (a, ao) = In A + An (/ (cc)) da (26) 


and ao is the initial conversion considered. 

Finally, this class of kinetic methods also includes 
the non-parameter kinetic (NPK) models (Sewry & 
Brown, 2002). These models assume that the reaction 
rate of a solid state reaction can be expressed as a 
product of two independent functions: one function 
depending on conversion and usually denoted as g(a) 
and the other depending on temperature and denoted 
as f(T). It should be noted that, in this case, the form 
of f(T) can be of a type other than the commonly used 
Arrhenius law type. This approximation is denoted as 
“single step” by Simon (2005), and has the mathemat¬ 
ical form expressed by the following equation 

P^ = g(a)f(T) (27) 

Several integral methods have already been stan¬ 
dardised for the experimental use of TGA and DSC 


Discussion 

Figs. 2a and 2b present the corresponding appar¬ 
ent activation energy and pre-exponential factor lit¬ 
erature values previously listed in Tables 1-3. Only 
two sets (A and E) of the kinetic parameters listed in 
Table 3 are shown. The results are coherent but very 
scattered, confirming the influence of the heat transfer 
and experimental conditions (type of cellulose, mass 
of sample, heating rate, pyrolysis temperature, etc.) 
on the reaction kinetics. Regarding the values of the 
activation energy, Varhegyi et al. (1989) state that the 
highest values are the most reliable. 

Fig. 2a shows that the highest values of the ki¬ 
netic parameters from the first order kinetic model 
(FI) are obtained when using microcrystalline cellu¬ 
lose (MCC). Other types of cellulose also afford the 
same results. These were not identified in the refer¬ 
ences chosen but it is possible they also were MCC. 
The same order of magnitude for these parameters 
is obtained using isoconversional methods, as is also 
the case for MCC. However, these values are situ¬ 
ated within the medium range of those obtained by 
using the FI model. With regard to the isoconver¬ 
sional methods, it should be noted that the Friedman 
method gives higher values for the activation energy 
than the Capela-Ribeiro method. 

Fig. 2b indicates that, in the case of the 3R model 
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applied to a sample of a-cellulose, the activation en¬ 
ergy required for the simultaneous formation of the 
three reaction products (char, tar, gases) is slightly 
higher than thatrequired for the formationof char: 
~ 200 kJ mol -1 vs ~ 170 kJ mol -1 . When comparing 
the values of E rl for reaction Eq. (8) from two models 
using two different types of cellulose, Fig. 2b indicates 
that higher values are obtained for the a-cellulose (3R 
model, ~ 200 kJ mol -1 ) than for the filter paper (BSZ 
model, ~ 150 kJ mol -1 ). 

It should also be noted that, for models BS, KB1, 
and DB, the values of the kinetic parameters are very 
close to each other, especially for reactions Eqs. (4), 
(10), and (12). Moreover, model AF exhibits slightly 
different values of these parameters when using dif¬ 
ferent types of cellulose. Hence, reactions Eqs. (17) 
and (18) give values of activation energy in the range 
of 130-150 kJ mol -1 when using a-cellulose and of 
approximately 170 kJ mol -1 when using filter pa¬ 
per. The corresponding values of the pre-exponential 
factor range are: from 1.94 x 10 11 min -1 up to 
1.06 x 10 13 min -1 for a-cellulose and approximately 
5.75 x 10 14 min -1 for filter paper. The reaction for 
the formation of intermediate from cellulose (reac¬ 
tion Eq. (4)) covers a wide range of values for the 
kinetic parameters. Thus, for the activation energy, 
these values are from approximately 130 kJ mol -1 up 
to 270 kJ mol -1 . The pre-exponential factor values are 
between 1.78 x 10 8 min -1 and 7.58 x 10 23 min -1 , re¬ 
spectively. The same behaviour can be observed for 
the reaction of tar formed from cellulose (reaction Eq. 
(2)): activation energy of 153-290 kJ mol -1 and pre¬ 
exponential factor of 1.58 x 10 15 -5.62 x 10 24 min -1 . 

Regarding the reaction of intermediate-forming 
(reaction Eq. (4)) that remains controversial, differ¬ 
ent behaviours can be observed when using differ¬ 
ent kinetic models. Hence, from the MKB model, 
the value of activation energy for this reaction is 
lower than for the reaction of tar-forming from cel¬ 
lulose (reaction Eq. (2)). At the same time, this 
value is the highest when using BS, DB, and DKM 
models. It is superior to the E a of the intermedi¬ 
ate decomposition to the pyrolysis products. This 
renders this reaction very slow in comparison with 
the others. Consequently, the intermediate-forming 
could be ignored in the kinetic models for this pro- 

Conclusions 

This paper presents a short review of the main ki¬ 
netic models published for cellulose pyrolysis. All of 
these models are to varying extents empirical, hence 
their degree of validity is generally limited to the par¬ 
ticular conditions used in their development (qual¬ 
ity of cellulose, heating rate, nature, flow-rate of car¬ 
rier gas, etc.). They are classified into two major 
types: “model-fitting” and “model-free” (isoconver- 
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Fig. 2. Compensation plots of several Arrhenius parameters from the literature (E a vs InA) for cellulose pyrolysis. FI (MCC (O), 
a-cellulose (x), filter paper (+), cellulose (^x)) and isoconversional (*) kinetic models (Tables 1 and 3) (a) and semi-global 
kinetic models KB1 (1975 ($), 1995 (☆)), BR (1998 (a)), BS (1979 (O), 1995 filter paper (•)), DKM (1987 (V)), 3R 
(a-cellulose: 1988 (A) 1995 (A)), MKB (a-cellulose: 1988 (<), 1995 (◄)), AF (1989: filter paper (>), a-cellulose (►), 1995 
filter paper (>)), DB (1994 (*)), KB (1997 (+)), BSZ (2001 filter paper (•)), F (2008 (x)) models (Table 2) (b). 

sional) methods. The first type uses a pre-defined perimental data on it while the second postpones the 
mathematical model for the process and fits the ex- choice of a model by first finding a variation in the ac- 
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tivation energy with the degree of conversion for the 
process studied. 

The kinetic models most frequently used for cel¬ 
lulose pyrolysis have been the FI and BS (Broido- 
Shafizadeh) models. The latter currently acts as a ref¬ 
erence for this process. When comparing the kinetic 
parameters values obtained from all these studies in 
the literature, a significant scattering could be ob¬ 
served but, nevertheless, a coherence. The differences 
are explained by the dissimilar experimental condi¬ 
tions and the lack of physical meaning of rate expres¬ 
sions and parameters, as observed by Galwey (2004). 
When using the FI model, the highest values of the 
kinetic parameters were obtained in the case of micro¬ 
crystalline cellulose (MCC). 

From a comparison of several values of these pa¬ 
rameters, when using semi-global kinetic models, it 
was concluded that the formation of an intermediate 
could be ignored due to the very high values obtained 
for the rate of this reaction. 

The author suggests that the FI kinetic model 
could be used for the design of a pyrolysis industrial 
plant, as a first approximation and for a general view 
of the process. For more accurate calculus, some iso- 
conversional method must be used. 

Acknowledgements. The author wishes to express her grat¬ 
itude to Mr. Yannick Soudais (Ecole des Mines d’Albi- 
Carmaux, France), Miss Maria-Octavia Lupgea (Centre Sci- 
entifique et Technique du Batiment - Grenoble, France) and 
Mr. Petru Budrugeac, for their help with the articles reviewed 
in this study, and Mr. Nicolae Sdrula (S.C. IPROCHIM S.A., 
Romania) and Mr. Gheorghe Soare (“Politehnica” University 
of Bucharest, Faculty of Applied Chemistry and Materials Sci¬ 
ence, Romania) for their valued comments and suggestions. 
The author would also like to thank the reviewers of the Chem¬ 
ical Papers Journal for their valuable commentaries and sug¬ 
gestions. 

References 

Aboulkas, A., & El Harfi, K. (2008). Study of the kinetics 
and mechanisms of thermal decomposition of Moroccan Tar- 
faya oil shale and its kerogen. Oil Shale, 25, 426-443. DOI: 
10.3176/oil.2008.4.04. 

Alves, S. S., & Figueiredo, J. L. (1989). Kinetics of cellulose 
pyrolysis modelled by three consecutive first-order reactions. 
Journal of Analytical and Applied Pyrolysis, 17, 37-46. DOI: 
10.1016/0165-2370(89)85004-1. 

AGEE-Stat (2011). Renewable energy sources 2010. Bonn, Ger¬ 
many: Federal Ministry for the Environment, Nature Conser¬ 
vation and Nuclear Safety. 

Agrawal, R. K. (1988). Kinetics of reactions involved in py¬ 
rolysis of cellulose. II. The modified Kilzer-Broido model. 
The Canadian Journal of Chemical Engineering, 66, 413- 
418. DOI: 10.1002/cjce.5450660310 
Banyasz, J. L., Li, S., Lyons-Hart, J. L., & Shafer, K. H. (2001). 
Cellulose pyrolysis: the kinetics of hydroxyacetaldehyde evo¬ 
lution. Journal of Analytical and Applied Pyrolysis, 57, 223- 
248. DOI: 10.1016/s0165-2370(00)00135-2. 

Barud, H. S., Ribeiro, C. A., Capela, J. M. V., Crespi, M. S., 
Ribeiro, S. J. L., & Messadeq, Y. (2011). Kinetic parame¬ 
ters for thermal decomposition of microcrystalline, vegetal, 
and bacterial cellulose. Journal of Thermal Analysis and 


Calorimetry, 105, 421-426. DOI: 10.1007/sl0973-010-1118- 
9. 

Blazek, J. (2005). Study of the reaction kinetics of the ther¬ 
mal degradation of polymer. Ph. D. thesis, Institut National 
Polytechnique de Toulouse, Toulouse, France. 

Broido, A., Javier-Son, A. C., Ouano, A. C., Barrell, E. 
M., II (1973). Molecular weight decrease in the early 
pyrolysis of crystalline and amorphous cellulose. Jour¬ 
nal of Applied Polymer Science, 17, 3627-3635. DOI: 
10.1002/app. 1973.070171207. 

Broido, A., & Nelson, M. A. (1975). Char yield on pyroly¬ 
sis of cellulose. Combustion and Flame, 24, 263-268. DOI: 
10.1016/0010-2180(75)90156-x. 

Broido, A. (1976). Kinetics of solid-phase cellulose pyrolysis. 
In F. Shafizadeh, K. V. Sarkanen, & D. A. Tillman (Eds.), 
Thermal uses and properties of carbohydrates and lignins 
(pp. 19-36). New York, NY, USA: Academic Press. DOI: 
10.1016/b978-0-12-637750-7.50006-6. 

Brown, M. E., Maciejewski, M., Vyazovkin, S., Nomen, R., 
Sempere, J., Burnham, A., Opfermann, J., Strey, R., An¬ 
derson, H. L., Kemmler, A., Keuleers, R., Janssens, J., 
Desseyn, H. O., Li, C. R., Tang, T. B., Roduit, B., Malek, J., 
& Mitsuhashi, T. (2000). Computational aspects of kinetic 
analysis. Part A: The ICTAC kinetics project-data, meth¬ 
ods and results. Thermochimica Acta, 355, 125-143. DOI: 
10.1016/s0040-6031(00)00443-3. 

Budrugeac, P., & Segal, E. (2003). Prediction of the isother¬ 
mal behavior of solid-gas systems from non-isothermal data. 
Journal of Thermal Analysis and Calorimetry, 72, 831—837. 
DOI: 10.1023/a:1025014114527. 

Cabrales, L., & Abidi, N. (2010). On the thermal degradation of 
cellulose in cotton fibers. Journal of Thermal Analysis and 
Calorimetry, 102, 485-491. DOI: 10.1007/sl0973-010-0911-9. 

Capart, R., Khezami, L., & Burnham, A. K. (2004). Assess¬ 
ment of various kinetic models for the pyrolysis of a micro- 
granular cellulose. Thermochimica Acta, 417, 79-89. DOI: 
10.1016/j.tca.2004.01.029. 

Chen, W. H., & Kuo, P. C. (2011). Isothermal torrefaction 
kinetics of hemicellulose, cellulose, lignin and xylan using 
thermogravimetric analysis. Energy, 36, 6451-6460. DOI: 
10.1016/j.energy.2011.09.022. 

Conesa, J. A., Caballero, J. A., Marcilla, A., & Font, R. (1995). 
Analysis of different kinetic models in the dynamic pyrol¬ 
ysis of cellulose. Thermochimica Acta, 254, 175-192. DOI: 
10.1016/0040-6031(94)02102-t. 

Di Blasi, C. (1993). Modeling and simulation of combus¬ 
tion processes of charring and non-charring solid fuels. 
Progress in Energy and Combustion Science, 19, 71-104. 
DOI: 10.1016/0360-1285(93)90022-7. 

Di Blasi, C. (1996). Heat transfer mechanisms and multi- 
step kinetics in the ablative pyrolysis of cellulose. Chemi¬ 
cal Engineering Science, 51, 2211—2220. DOI: 10.1016/0009- 
2509(96)00078-4. 

Di Blasi, C. (1998). Comparison of semi-global mechanisms for 
primary pyrolysis of lignocellulosic fuels. Journal of Analyt¬ 
ical and Applied Pyrolysis, 47, 43—64. DOI: 10.1016/s0165- 
2370(98)00079-5. 

Dickinson, C. F., & Heal, G. R. (2009a). A review of the IC¬ 
TAC Kinetics Project, 2000. Part 1. Isothermal results. Ther¬ 
mochimica Acta, 494, 1-14. DOI: 10.1016/j.tca.2009.05.003. 

Dickinson, C. F., & Heal, G. R. (2009b). A review of the ICTAC 
kinetics project, 2000. Part 2. Non-isothermal results. Ther¬ 
mochimica Acta, 494, 15-25. DOI: 10.1016/j.tca.2009.05.009. 

Diebold, J. P. (1994). A unified, global model for the pyrol¬ 
ysis of cellulose. Biomass and Bioenergy, 7, 75-85. DOI: 
10.1016/0961-9534(94)00039-v. 






C. Serbariescu Chemical Papers 68 (7) 847-860 (2014) 859 


Evans, R. J., & Milne, T. A. (1987). Molecular characterization 
of the pyrolysis of biomass. 1. Fundamentals. Energy & Fuels, 
1, 123-137. DOI: 10.1021/ef00002a001. 

Fisher, T., Hajaligol, M., Waymack, B., & Kellogg, D. (2002). 
Pyrolysis behavior and kinetics of biomass derived materials. 
Journal of Analytical and Applied Pyrolysis, 62, 331—349. 
DOI: 10.1016/s0165-2370(01)00129-2. 

Font, R., & Garcia, A. N. (1995). Application of the transi¬ 
tion state theory to the pyrolysis of biomass and tars. Jour¬ 
nal of Analytical and Applied Pyrolysis, 35, 249-258. DOI: 
10.1016/0165-2370(95)00916-8. 

Galwey, A. K. (2004). Is the science of thermal analysis kinet¬ 
ics based on solid foundations? A literature appraisal. Ther- 
mochimica Acta, 413, 139-183. DOI: 10.1016/j.tca.2003.10. 
013. 

Garcia-Perez, M. (2008). The formation of polyaromatic hydro¬ 
carbons and dioxins during pyrolysis: A review of the litera¬ 
ture with descriptions of biomass composition, fast pyrolysis 
technologies and thermochemical reactions. Pullman, WA, 
USA: Washington State University. (WSUEEP08-010) 

Gavillon, R. (2007). Preparation et caracterisation des materi- 
aux cellulosiques ultra poreux. Ph. D. thesis, iScole des Mines 
de Paris, Paris, France, (in French) 

Grpnli, M., Antal, M. J., Jr., & Varhegyi, G. (1999). A round- 
robin study of cellulose pyrolysis kinetics by thermogravime¬ 
try. Industrial & Engineering Chemistry Research, 38, 2238- 
2244. DOI: 10.1021/ie980601n. 

Hopkins, M. W., DeJenga, C., & Antal, M. J., Jr. (1984). The 
flash pyrolysis of cellulosic materials using concentrated vis¬ 
ible light. Solar Energy, 32, 547-551. DOI: 10.1016/0038- 
092x(84)90269-x. 

Hu, S., Jess, A., & Xu, M. H. (2007). Kinetic study of Chinese 
biomass slow pyrolysis: Comparison of different kinetic mod¬ 
els. Fuel, 86, 2778-2788. DOI: 10.1016/j.fuel.2007.02.031. 

Huang, Y. F., Kuan, W. H., Chiueh, P. T., & Lo, S. L. (2011). 
A sequential method to analyze the kinetics of biomass 
pyrolysis. Bioresource Technology, 102, 9241-9246. DOI: 
10.1016/j.biortech.2011.07.015. 

Kim, S. D., & Eom, Y. J. (2006). Estimation of kinetic triplet 
of cellulose pyrolysis reaction from isothermal kinetic results. 
Korean Journal of Chemical Engineering, 23(3), 409-414. 
DOI: 10.1007/bf02706742. 

Kilzer, F. J., & Broido, A. (1965). Speculations on the nature 
of cellulose pyrolysis. Pyrodynamics, 2, 151-163. 

Lede, J. (2012). Cellulose pyrolysis kinetics: An historical review 
on the existence and role of intermediate active cellulose. 
Journal of Analytical and Applied Pyrolysis, 94, 17-32. DOI: 
10.1016/j.jaap.2011.12.019. 

Lewellen, P. C., Peters, W. A., & Howard, J. B. (1977). Cellu¬ 
lose pyrolysis kinetics and char formation mechanism. Sym¬ 
posium (International) on Combustion, 16, 1471-1480. DOI: 
10.1016/s0082-0784(77)80429-3. 

Li, C. R., & Tang, T. B. (1997). Dynamic thermal analy¬ 
sis of solid-state reactions. The ultimate method for data 
analysis? Journal of Thermal Analysis, 49, 1243—1248. DOI: 
10.1007/bf01983680. 

Li, C. R., & Tang, T. B. (1999). Isoconversional method for ki¬ 
netic analysis of solid-state reactions from dynamics thermo- 
analytical data. Journal of Materials Science, 34, 3467-3470. 
DOI: 10.1023/a: 1004605820783. 

Liao, Y. F., Wang, S. R., & Ma, X. Q. (2004). Study of reac¬ 
tion mechanisms in cellulose pyrolysis. Preprints of Papers- 
American Chemical Society, Division of Fuel Chemistry, 49, 
407-411. 

Liu, Q., Wang, S. R., Wang, K. G., Guo, X. J., Luo, Z. Y., 
& Cen, K. F. (2008). Mechanism of formation and con¬ 
sequent evolution of active cellulose during cellulose py¬ 


rolysis. Acta Physico-Chimica Sinica, 24, 1957-1963. DOI: 
10.1016/sl872-1508(08)60078-9. 

Mamleev, V., Bourbigot, S., & Yvon, J. (2007). Kinetic analysis 
of the thermal decomposition of cellulose: The change of the 
rate limitation. Journal of Analytical and Applied Pyrolysis, 
80, 141-150. DOI: 10.1016/j.jaap.2007.01.012. 

Mamleev, V., Bourbigot, S., Le Bras, M., & Yvon, J. (2009). 
The facts and hypotheses relating to phenomenological model 
of cellulose pyrolysis: Interdependence of the steps. Jour¬ 
nal of Analytical and Applied Pyrolysis, 84, 1-17. DOI: 
10.1016/j.jaap.2008.10.014. 

Marra, F. (2009). Numerical analysis for kinetics and yield of 
wood biomass pyrolysis. In N. Mastorakis, & J. Sakellaris 
(Eds.), Advances in numerical methods (chapter 11, pp. 127- 
136). Heidelberg, Germany: Springer. DOI: 10.1007/978-0- 
387-76483-211. 

Miller, R. S., & Bellan, J. (1997). A generalized biomass pyroly¬ 
sis model based on superimposed cellulose, hemicellulose and 
lignin kinetics. Combustion Science and Technology, 126, 97- 
137. DOI: 10.1080/00102209708935670. 

Ranzi, E., Cuoci, A., Faravelli, T., Frassoldati, A., Migliavacca, 
G., Pierucci, S., & Sommariva, S. (2008). Chemical kinetics 
of biomass pyrolysis. Energy & Fuels, 22, 4292-4300. DOI: 
10.1021 / ef800551t. 

Reed, T. B., & Cowdery, C. D. (1987). Heat flux requirements 
for fast pyrolysis and a new method for generating biomass 
vapor. In 193rd National Meeting of the American Chemi¬ 
cal Society, April 5-10, 1987. Denver, CO, USA: American 
Chemical Society Division of Petroleum Chemistry. 

Sanchez-Jimenez, P. E., Perez-Maqueda, L. A., Perejon, A., 
Pascual-Cosp, J., Benitez-Guerrero, M., & Criado, J. M. 
(2011). An improved model for the kinetic description of the 
thermal degradation of cellulose. Cellulose, 18, 1487-1498. 
DOI: 10.1007/sl0570-011-9602-3. 

Sbirrazzuoli, N., Vincent, L., Mija, A., & Guigo, N. (2009). 
Integral, differential and advanced isoconversional methods. 
Complex mechanisms and isothermal predicted conversion¬ 
time curves. Chemometrics and Intelligent Laboratory Sys¬ 
tems, 96, 219-226. DOI: 10.1016/j.chemolab.2009.02.002. 

Scott, D. S., Piskorz, J., & Radlein, D. (1989). Thermal con¬ 
version of biomass to liquids by the Waterloo fast pyrolysis 
process. In E. Mattucci, G. Grassi, & W. Palz (Eds.), Pro¬ 
ceedings of Pyrolysis as a Basic Technology for Large Agro- 
Energy Projects, October 15 Hi. 1987 (pp. 115-124). Brus¬ 
sels, Belgium: Office for Official Publications of the European 
Communities. 

Sewry, J. D., & Brown, M. E. (2002). “Model-free” ki¬ 
netic analysis? Thermochimica Acta, 390, 217-225. DOI: 
10.1016/s0040-6031(02)00083-7. 

Shafizadeh, F. (1968). Pyrolysis and combustion of cellulosic 
materials. Advances in Carbohydrate Chemistry, 23, 419- 
474. DOI: 10.1016/s0096-5332(08)60173-3. 

Shafizadeh, F., & Bradbury, A. G. W. (1979). Thermal degra¬ 
dation of cellulose in air and nitrogen at low temperatures. 
Journal of Applied Polymer Science, 23, 1431-1442. DOI: 
10.1002/app. 1979.070230513. 

Shafizadeh, F. (1982). Introduction to pyrolysis of biomass. 
Journal of Analytical and Applied Pyrolysis, 3, 283-305. 
DOI: 10.1016/0165-2370(82)80017-x. 

Simon, P. (2005). Considerations on the single-step kinetics ap¬ 
proximation. Journal of Thermal Analysis and Calorimetry, 
82, 651-657. DOI: 10.1007/sl0973-005-0945-6. 

Sonobe, T., & Worasuwannarak, N. (2008). Kinetic analyses 
of biomass pyrolysis using the distributed activation energy 
model. Fuel, 87, 414-421. DOI: 10.1016/j.fuel.2007.05.004. 

Stamm, A. J. (1956). Thermal degradation of wood and cellu¬ 
lose. Industrial & Engineering Chemistry, 48, 413-417. DOI: 
10.1021/ie51398a022. 




C. 5erbar 


j/Chemical Papers 68 (7) 847-860 (2014) 


Tihay, V., Boulnois, C., & Gillard, P. (2011). Influence of oxygen 
concentration on the kinetics of cellulose wadding degrada¬ 
tion. Thermochimica Acta, 525, 16-24. DOI: 10.1016/j.tca. 
2011.07.016. 

Varhegyi, G., Antal, M. J., Jr., Szekely, T., Szabo, P. (1989). 
Kinetics of the thermal decomposition of cellulose, hemicel- 
lulose, and sugarcane bagasse. Energy & Fuels, 3, 329-335. 
DOI: 10.1021/ef00015a012. 

Varhegyi, G., Jakab, E., & Antal, M. J., Jr. (1994). Is the 
Broido-Shafizadeh model for cellulose pyrolysis true? Energy 
& Fuel, 8, 1345-1352. DOI: 10.1021/ef00048a025. 

Varhegyi, G., Antal, M. J., Jr., Jakab, E., & Szabo, P. (1997). 
Kinetic modeling of biomass pyrolysis. Journal of Analyti¬ 
cal and Applied Pyrolysis, 42, 73-87. DOI: 10.1016/s0165- 
2370(96)00971-0. 

Volker, S., & Rieckmann, Th. (2002). Thermogravimetric in¬ 
vestigation of cellulose pyrolysis - impact of initial and final 
mass on kinetic results. Journal of Analytical and Applied 
Pyrolysis, 62, 165-177. DOI: 10.1016/s0165-2370(01)00113- 
9. 

Vyazovkin, S. (1996). A unified approach to kinetic pro¬ 
cessing of nonisothermal data. International Journal of 
Chemical Kinetics, 28(2), 95-101. DOI: 10.1002/(SICI)1097- 
4601(1996)28:2<95::AID-KIN4>3.0.CO;2-G. 


Vyazovkin, S., & Dollimore, D. (1996). Linear and nonlin¬ 
ear procedures in isoconversional computations of the ac¬ 
tivation energy of nonisothermal reactions in solids. Jour¬ 
nal of Chemical Information and Modeling, 36, 42-45. DOI: 
10.1021/ci950062m. 

Vyazovkin, S., & Wight, C. A. (1997). Kinetics in solids. 
Annual Review of Physical Chemistry, 48, 125-149. DOI: 
10.1146/annurev.physchem.48.1.125. 

White, J. E., Catallo, W. J., & Legendre, B. L. (2011). Biomass 
pyrolysis kinetics: A comparative critical review with relevant 
agricultural case studies. Journal of Analytical and Applied 
Pyrolysis, 91, 1-33. DOI: 10.1016/j.jaap.2011.01.004. 

Zhu, G. Y., Zhu, X., Xiao, Z. B., & Yi, F. P. (2012). 
Study of cellulose pyrolysis using an in situ visualiza¬ 
tion technique and thermogravimetric analyzer. Journal 
of Analytical and Applied Pyrolysis, 94, 126-130. DOI: 
10.1016/j.jaap.2011.11.016. 



